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The surface spin states for bismuth thin films are investigated using an sp^ tight-binding 
model. The model explains most experimental observations using angle-resolved photoemission 
spectroscopy, including the Fermi surface, the band structure with Rashba spin splitting, and the 
quantum confinement in the energy band gap of the surface states. A large out-of-plane spin com¬ 
ponent also appears. The surface states penetrate inside the film to within approximately a few 
bilayers near the Brillouin-zone center, whereas they reach the center of the film near the Brillouin- 
zone boundary. 


Introduction. — The spin-orbit interaction (SOI) in¬ 
duces spin splitting in the absence of an external mag¬ 
netic field on a two-dimensional (2D) system, i.e., Rashba 
spin splitting^, which has been an indispensable element 
of spintronic physics and devices^. The Rashba effect is 
expected on crystal surfaces due to their inversion asym¬ 
metry. For example, Rashba spin splitting has been ob¬ 
served on the Au(lll) surface^”^. Bismuth (Bi) is a 
group V semimetal with a large SOI due to the heavy 
mass of the Bi atom; therefore, the surface of Bi crystals 
is an ideal system to observe a strong Rashba effect^. 

Angle-resolved photoemission spectroscopy (APRES) 
experiments have been reported for the Bi surface accom¬ 
panied with first-principles band calculations^^^^. The 
surface states have a hexagonal electron pocket around 
the F point and six-fold hole pockets 7,9,io,i6,i9 
principles band calculations showed that these two sur¬ 
face states are spin-split bands and this Rashba split¬ 
ting has been confirmed experimentally In addi¬ 

tion, the surface spin orientation has been elucidated, 
and in particular, a giant out-of-plane spin polarization 
was reported The band structure is dependent on 
the film thickness because of the quantum confinement 
effect^^’^^’^^. 

In addition to the ARPES experiments, many inter¬ 
esting features have been studied in the electronic trans¬ 
port properties of Bi nanostructures. The conductivity 
of Bi films has been measured , and was determined 
to be dependent on both the surface and bulk states, 
and the coupling between them has a major influence on 
the conductivity Quantum confinement effects can sig¬ 
nificantly enhance thermoelectric properties in quantum 
well and quantum wire structures'^. Bismuth is thus 
a prime candidate to achieve high performance thermo¬ 
electric conversion in its nanostructures^^’^^. To under¬ 
stand the transport properties in Bi quantum confine¬ 
ment structures, it is necessary to simultaneously deter¬ 
mine the electronic properties of both the surface and the 
bulk states. 

Although first-principles band calculations have al¬ 
ready revealed the Eermi surface and the energy band 
structure 8 , 10 , 12 , 26 ^ systematic analysis for compari¬ 
son with the reported ARPES experimental results has 
been conducted to date. Here, we approach this issue us¬ 


ing an sp^ tight-binding model that reproduces the band 
structure of bulk Bi proposed by Liu and Allen^^. This 
model has been applied to discuss the topological and 
non-topological phases of the surface states of pure Bi 
and Sb^^, and Bii-a^Sb^^^^, as well as two-dimensional 
Bi^^. Extra surface hopping terms are added that 
were originally proposed to explain the Au(III) surface 
states. This model will enable confirmation of whether 
the ARPES results originate from the surface effect. In 
addition, it is straightforward to see the effects of quan¬ 
tum confinement because the film thickness can be easily 
changed and the electronic states both inside the film 
and at the surface can be analyzed, which is important 
to investigate the electronic transport properties. We can 
thus give a systematic survey of the ARPES experimental 
results by taking advantage of these points. 

Model Hamiltonian. — Bismuth has a rhombohedral 
Bravais lattice with two atoms per unit cell, forming a 
bilayer (BL) structure, as shown in Eig. 1(a). A Bi thin 
film is obtained by stacking the BLs along the (III) di¬ 
rection, such as the z-axis depicted in Eig. 1(b). The 
surface is thus parallel to the xy plane. The uppermost 
and lowermost BLs are in contact with a vacuum. 

We first construct a model Hamiltonian for the Bi thin 
film. Eor this purpose, the sp^ tight-binding model pro¬ 
posed for the bulk Bi crystal^^ is adapted to the Bi thin 
film. There are 5 -, Px-, Py-^ and p^-orbitals with spin 
index a on each atom. The hopping terms among the 
atomic orbitals are decomposed into inter- and intra-BL 
hopping terms. The inter-BL hopping term 7721-2 con¬ 
sists of the nearest-neighbor hopping term in the bulk Bi 
Hamiltonian, whereas the intra-BL hopping term consists 
of two parts, Hu and 77i2-i, with the third and second 
nearest-neighbor hopping terms in the bulk model, re¬ 
spectively. The Eermi energy is set to zero. 

There is a surface potential gradient on the surface 
BL along the z-axis between the surface Bi atoms and 
the vacuum. The surface Rashba effect is induced by 
the contribution of this potential gradient^^’^^. In terms 
of the sp^ tight-binding model, this is described by the 
following spin independent hopping terms between the 
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FIG. 1. (Color online) (a) Bismuth crystal structure, (b) Schematic view of Bi thin film, (c, d) Surface BL density of states 
(DOS) at the Fermi energy /?(!, 0,/c^) of a 16 BL film with (c) 7pp/sp = 0, and (d) with 7 sp = 0.45 and 7 pp = —0.27. (e, 
f) Energy band structure along the F-M line {kx = 0). The DOS p{i, E,0,ky) is plotted for the (e) middle BL {i = 8) and 
(f) the surface BL (z = 1). The white lines represent the eigenvalues of H. (g) Spin-resolved band structure for (f), where 
s = Sx{l, E, 0, ky). Note that Sy/z(l, E, 0, ky) = 0. 


nearest-neighbors sites, Ri and Rj 


represented by the following matrix form: 



7pp ^ij 


iPx,Pz) or {Pz,Px), 


( -f^sll 

Hi2-2 

H21-2 

H12-1 

\ 

^a/3 — ^ 

7 pp sin Oij 

,(q:,/3) = 

(jpy,Pz) or iPz,Py), (1) 

H = 



H21-1 

Hii 



^7sp 

,(q:,/3) = 

{s,pz) or {pz,s), 


\ 




Hii iZi 2 -i 
H21-1 ff'n/ 


where Oij is the azimuthal angle between Ri — Rj and the 
x-axis, and 7 pp and 7sp are the hopping matrix elements 
of the Hamiltonian. Note that those hopping terms are 
zero in the bulk Bi crystal model because of the inver¬ 
sion symmetry. It is assumed that the surface hopping 
terms (1) appear only on the uppermost atomic layer, 
the first atomic layer of the uppermost BL, and the low¬ 
ermost atomic layer with —7sp/pp, because the surface 
field points in the opposite direction at the lowermost 
layer. The values of 7sp and 7 pp remain to be determined 
at this stage. 


( 2 ) 

where iLgii is the Hamiltonian for the uppermost atomic 
layer, which includes the surface hopping terms (1) in 
addition to iLn, while i^ Hamiltonian that in¬ 
cludes the surface hopping terms (1) with —7pp/sp. The 
size of the matrix is thus 16n x 16n when the number of 
the BLs is n. The Hamiltonian (2) is a function of the 
wave vectors kx and ky: H = H{kx^ky). 

Calculation of DOS .— The DOS and the band struc¬ 
ture of the thin film is obtained from the retarded Green’s 
function matrix G{E,kx^ky) with energy E defined by 


The total Hamiltonian of the thin film H is therefore 


G{E,kx,ky) = [E + i6 - H{kx,ky)] ^ 


(3) 
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with (5 = 1.0 X 10 ^ in the numerical calculations. The 
DOS in the ith BL is defined by 

p{i,E,kx,ky) = --TTlmG{E,kx,ky), (4) 

TT 

where Tr represents the trace over the orbitals and the 
spin only on the ith BL {i = 1 for the uppermost BL). In 
a similar way, the spin-resolved DOS Sa{i, E^ kx^ ky) {a = 
X, y, z) is given by 

^— Tr ImS qjCt(L/, ), (5) 

TT 

where Sq, is the Pauli spin matrix that acts on the four 
orbital states. The eigenvalues of H are also calculated 
to show the entire band structure of the film. 

Parameter fitting. — In the following, 7sp and 7pp are 
treated as fitting parameters. To fix these values, we 
use a phenomenological approach: We first calculate the 
DOS on the surface BL and the band structure for various 
values of 7 sp/pp and then compare them with the ARPES 
experimental results 7,8,io,ii,i6,i9 ^^est selec¬ 

tion. This scheme was successful and led to 7sp = 0.45, 
and 7 pp = —0.27. The numerical results appear similar 
near these values. Note that these values are the same 
order of magnitude as the hopping matrix elements be¬ 
tween the second and third nearest neighbors given in 
Ref.^A 

The presence of the surface terms (1) is essential to ex¬ 
plain the observed Fermi surface. Figures 1(c) and (d), 
show the DOS on the surface BL at the Fermi energy 
p{l,0,kx^ky) for the 16 BL thin film without and with 
the surface hopping terms, respectively. In both cases, 
a hexagonal electron pocket appears around the f point 
designated by Si. Qualitative differences arise outside 
of Si] with the surface hopping term, there are six hole 
lobes and six extra electron lobes, designated by S' 2 , and 
S's, respectively, while S 2 is missing without the surface 
hopping term. The ARPES experiments show the pres¬ 
ence of S 2 , which confirms that the surface terms (1) play 
a central role in the formation of the Fermi surface. 

Band strueture. — Next, we discuss the energy band 
structure along the F-M line {kx = 0). Figures 1(e) and 
(f) show p{i, E^O^ky) for the middle {i = 8) and sur¬ 
face (i = 1) BLs, respectively. The eigenvalues of H are 
also shown as white lines for comparison. On the mid¬ 
dle BL, the plot covers most of the eigenvalues, while on 
the surface BL, the plot appears only in a small frac¬ 
tion of the eigenvalue curves and mostly on two curves 
near the Fermi energy. The upper curve forms the Si 
and S's structures, whereas the lower curve forms the S 2 
structure. 

The spin-resolved band structure illustrates the dis¬ 
tinctive features of the surface states, as shown in 
Fig 1(g). The spin splitting appears near the F point, 
which is similar to Rashba spin splitting, and it dimin¬ 
ishes near M. This is consistent with the experimen¬ 
tal results and the first-principles band calculations^^’^^. 
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FIG. 2. (Color online) (a) In-plane surface spin magnitude 
|S|||. The green arrows indicate the direction of s,, for repre¬ 
sentative points, (b) Out-of-plane surface spin s^. The region 
surrounded by a dotted line corresponds to that shown in (a). 
Inset: Sz around the F point without the surface hopping 
terms and with the same scale as the main figure. Additional 
information on the surface spin states is given in the Supple¬ 
mental MateriaP^. 


Thus, the surface states on the Bi film are well described 
by the phenomenological tight-binding model. 

Surfaee spin states. — Next, we discuss the surface spin 
texture at the Fermi energy; s^ = Sa{l,0,kx,ky). Fig¬ 
ure 2(a) shows the in-plane spin = (sx^Sy) distribu¬ 
tion. On Si, S|| lies along the pocket structure, while on 
S 2 ^ the direction of s is opposite to that on Si. The in¬ 
plane spin rotations on Si and S 2 are broadly similar to 
those by the Rashba SOL In addition, |S||| along the lobe 
on S 2 is almost constant. These observations are con¬ 
sistent with previous experimental results^’^^’^^’^^. How¬ 
ever, the asymmetry of the |S||| along the ky axis on S 2 ^^ 
is not observed in the present model. Instead of this 
asymmetry, |S||| on Si oscillates every 60°. In addition, s 
is not always perpendicular to k on S 2 ^ which comes from 
the non-parabolic band structure. These indicate that 
the spin structure is not described by a simple Rashba 
SOI model. 

The deviation from the simple Rashba model is clar¬ 
ified by the out-of-plane spin Sz, as shown in Fig. 2(b). 
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FIG. 3. (Color online) (a) Surface BL DOS at the Fermi energy p(l, 0, fc, of a 40 BL film with two lines and three points, 
(b, c) Energy band structure for various numbers of BLs. The surface BL DOS p{l^E^kx^ky) is plotted along the lines (I) 
{ky = 0.22) and (II) {ky = 0.5), respectively, (d-i) The layer-resolved DOS at the Fermi energy, p{n) and Sx(n), at ky = 0.08, 
0.16, and 0.65, the points (A), (B), and (C) in (a), respectively, for (d-f) 16 BL, and (g-i) 100 BL. In (g) and (h), p(n) and 
Sxiji) are shown only near the surface BLs. 


There is a relatively large Sz over the SiSs structures, 
where the maximum of \sz\ is approximately 25% of the 
maximum of |S|||. Furthermore, Sz changes its sign every 
60°. These results are consistent with recent experimen¬ 
tal results^^, although larger values of \sz\ are observed 
experimentally. In addition, the fine structure of Sz is 
clarified, where the sign of Sz also changes from Si to S 3 
in the same manner as s^^. 

We further discuss the presence of the giant Sz com¬ 
ponent. The first-principles calculations show similar re¬ 
sults for Sz with the topological phase of the Bii_a,Sba; 
crystal^"^, where Sz is very small and around 1% of |S|||, 
which indicates the spin lies on the two-dimensional sur¬ 


face. Pure Bi {x = 0) is in the trivial phase^^’^^’^^; there¬ 
fore, large values of Sz may be direct evidence for clarifi¬ 
cation of the difference between the trivial and topologi¬ 
cal phases, besides the number of Fermi surface crossings 
from the zone center to the boundary. To support this 
point within the proposed model, Sz without the surface 
hopping term near the F point is shown in the inset of 
Fig. 2(b). Similar results are obtained both with and 
without the surface hopping terms. This indicates that 
the origin of Sz is not from the surface effect, but from 
the bulk hopping terms and the atomic SOI of Bi itself, 
which determines the bulk band structure. Thus, the 
origin of Sz for pure Bi is associated with the bulk band 
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structure, which leads to a trivial phase. 

BL number dependence .— Figure 3(a) shows 
p(l^{)^kx^ky) for a 40 BL film. Compared with 
that for the 16 BL film, both Si and S 2 structures are 
unchanged, while the structure is prolonged towards 
the M point, which is consistent with the experimental 
results^^’^^. To examine this difference in detail, we 
discuss the BL number dependence along the two lines 
and three points shown in Fig. 3(a). 

Figure 3(b) shows p{l^E^kx^ky) in S 2 for various num¬ 
bers of BLs along the line (I) shown in Fig. 3(a). The 
two surface states near the Fermi energy separated by a 
band gap are not affected by changing the BL number. 
However, the energy levels away from the Fermi level 
are under the strong influence of the BL number, which 
indicates the quantum confinement in the thin film. Fig¬ 
ure 3(c) shows p{l, E^kx^ky) in Ss along the line (II) 
shown in Fig. 3(a). Although the band structure near the 
Fermi level shows a linear dispersion similar to that for 
line (I), the band gap clearly decreases as the BL number 
increases. A similar observation is obtained experimen¬ 
tally^^. Hence, the surface states are under the strong 
influence of quantum confinement on Ss, while they are 
not on 5'2- 

Finally we discuss the surface state penetration in¬ 
side the thin film. For this purpose, the layer-resolved 
DOS, p{n) = p{n, 0, 0, ky) and Sx{n) = Sx{n, 0, 0, ky), are 
shown in Figs. 3(d-f) and (g-i) for the 16 BL and 100 
BL films, respectively, at the three points indicated in 
Fig. 3(a). All the figures show that the spin on the up¬ 


permost and lowermost BLs are in opposite directions, as 
expected; The spin changes its sign at the middle of the 
film. The surface states on Si penetrate only a few BLs, 
and a similar result is obtained for the surface states on 
S 2 with a slightly longer penetration length. The pene¬ 
tration length is unchanged by the film thickness, which 
confirms they are genuine surface states. On other hand, 
at S's, p{n) and Sx{n) decay over 20 BLs, and p(n) is 
finite even at the middle of the film. Thus, the states 
are no longer simple “surface” states and are under the 
influence of the quantum confinement inside the film. 

Conclusions .— We have shown that an sp^ tight- 
binding model with surface hopping terms can explain 
most of the experimental ARPES observations for bis¬ 
muth thin films, including the Fermi surface, the spin- 
resolved band structure with Rashba spin splitting, and 
the quantum confinement effect in the energy band struc¬ 
ture. The model also explains the large out-of-plane spin 
observed, which originates from the intrinsic Bi crystal 
structure rather than the surface effect. We have also 
clearly shown that the surface states penetrate inside 
the film to within approximately a few BLs near the 
Brillouin-zone center, whereas they reach the center of 
the film near the Brillouin-zone boundary. 
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Supplemental Material 

Tight-binding theory of surface spin states on bismuth thin films 

This supplemental material contains 

• Magnitude and direction of the in-plane spin, 

• Surface spin density of states (DOS) at the Fermi energy. 


Magnitude and direction of the in-plane spin. — Figure Sl(a) shows the detailed magnitude and direction of s^^ and 
Fig. Sl(b) shows a plot of |S||| for reference. The green arrows in Fig. 2(a) of the main text represent the directions 
of S|| shown in Fig. SI (a) as representative points. 

Surface spin DOS at the Fermi energy. — Figure S2 shows the surface spin DOS at the Fermi energy for all spin 
components. 
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FIG. SI. (a) Magnitude and direction of s,, as functions of the wave vectors kx and ky. The length of the arrow at each point 
is proportional to |S|||. (b) |S||| for reference. The same as Fig. 2(a) in the main text without the magnitude scale and green 
arrows. 
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FIG. S2. Surface spin DOS at the Fermi energy: Sa = Sa(l, 0, fc,/c^) with a = x,y,z. 










